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ABSTRACT 

We show that a star orbiting close enough to an adiabatically grown supermassive black hole (SMBH) can 
capture weakly interacting massive particles (WIMPs) at an extremely high rate. The stellar luminosity due to 
annihilation of captured WIMPs in the stellar core may be comparable to or even exceed the luminosity of the 
star due to thermonuclear burning. The model thus predicts the existence of unusual stars, essentially WIMP 
burners, in the vicinity of a SMBH. We find that the most efficient WIMP burners are stars with degenerate 
electron cores, e.g. white dwarfs (WDs); such WDs may have a very high surface temperature. If found, such 
stars would provide evidence for the existence of particle dark matter and can possibly be used to establish its 
density profile. On the other hand, the lack of such unusual stars may provide constraints on the WIMP density 
near the SMBH, as well as the WIMP-nucleus scattering and pair annihilation cross-sections. 
Subject headings: black hole physics — elementary particles — radiation mechanisms: non-thermal — stars: 
general — stars: evolution — dark matter 



1. INTRODUCTION 

The nature of the non-baryonic dark matter, which dom- 
inates the visible matter by about 4:1, is perhaps the most 
interesting experimental challenge for contemporary particle 
astrophysics. A hint for a solution has been found in particle 
physics where the WIMPs arise naturall y in supersymmetri c 
extensions of the Standard Model (e.g., Habe r & Kandll985h . 
among other possibilities. The WIMP is typically defined 
as a stable, electrically neutral, massive particle. Assuming 
that non-baryonic dark matter is dominated by WIMPs, the 
pair ann ihilation cross-section is related to the observed relic 
density ( Bergstrom 2000; lJungman. Kamion kowski & Griesj 
[1996). A pair of WIMPs can annihilate producing ordinary 
particles and 7-rays. 

WIMPs are expected to form high density clumps according 
to N-body simul ations of te st partic les with only gravitational 
intera ctions (Moor e et al.l 119991 : iNavarro. Frenk. & White! 
Il997h . The highest density "free space" dark mat- 
ter regions occur for dark matter particles captured 
within the gravitational pot ential of adiabatically grown 
SMBHs (He rtone & Merritt 200l iGnedin & Primackll2004 
iGondolo & Silki 1 999) . Higher dark matter densities are pos- 
sible for dark matter particles captured inside of stars or plan- 
ets. Any star close enough to a SMBH can capture a large 
number of WIMPs during a short period of time. Annihilation 
of captured WIMPs may lead to considerable energy release 
in stellar cores thus affecting the evolution and appearance of 
such stars. 

Such an idea has been fir st propos ed by Salati & Silk (il989h 
and further developed by iBouquet & Salati (1989.) who ap- 
plied it to main-sequence stars. The model led to the conclu- 
sion of suppression of stellar core convection, thus predicting 
a concentration of stars in the Galactic Center masquerading 
as cold red giants. 

An order-of-magnitude estimate of the WIMP capture 
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rates for stars of various masses and evolution stages 
dMoskalenko & Wail 12006!) lead us to the conclusion that 
WDs, fully burned stars without their own energy supply, are 
the most promising candidates to look for In this paper we 
calculate the WIMP capture by WDs located in a high den- 
sity dark matter region, and discuss their observational fea- 
tures. We use current limits on WIMP-nucleus interaction and 
WIMP annihilation cross sections, as well as recent estimates 
of WIMP energy density near an adiabatically grown SMBH. 

2. WIMP ACCUMULATION IN STARS 

In a steady state the WIMP capture rate C is balanced by 
the annihilation rate (j3riest & Secke L1987.) 



where 



C = ANt 



X 
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{(Jav) is the velocity averaged WIMP pair annihilation cross- 
section, the effective radius 
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is determined by matching the star core temperature Tc with 
the gravitational potential energy (assuming thermal equilib- 
rium), c is the speed of light, G is the gravitational constant, 
Pc is the star core density, and is the WIMP mass. The 
total number of WIMPs captured by a star is 

= CTeq tanh(T*/re,), (4) 

where is the star's age, and the equilibrium time scale is 
given by 

= (5) 

Limits from direct detection of dark matter on the WIMP- 
nucleon cross-section imply that only a fraction of the WIMPs 
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Fig. 1 . — WIMP capture rate by Oxygen WDs (An = 16) vs. velocity dispersion for 
envelope and T = 100, 000 K (Panel et al. 2000). Right panel corresponds to the zero-temperature approximation (Hamada & Salpetei 1961). The solid hnes 



pmax jjjg [gj!( panel shows calculations for WDs without hydrogen 
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show the capture rate calculated using the modified cross section (Eq. (9)), the dashed lines are calculated for CTq = (tq. The labels show WD mass in A/q units. 



crossing the star will scatter and be captured. The capture 
rate for a Maxwellian WIMP velocity distribution (in the ob- 
server's frame) by a star moving wi th an arbitrar y velocity w» 
relative to the observer is given bv (lGouldll987h : 



C = 47r 
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dr r 



where 

dC{r) 



dV 



1/2 



aoA, 



4 P* 



dV ' 



' Mn my 



(6) 



(7) 



1 



- 2 ) [xi-V,v)-xiA-,A+)] 



.-At 



-A-e ^+ 



rje 



A' = 



3i;2(r)^ 



A±=A±r], 



?7 = 



Xia.b): 



dye' 



[erf(6) -erf(a)], 



is the ambient WIMP energy density, An is the atomic 
number of the star's nuclei, M„ is the nucleus mass, v is the 
WIMP velocity dispersion, and /i — m^lMn, /i_ = (/i — 
l)/2. The escape velocity at a given radius r inside of a star 
is given by 



v{r) 



2G / dV 



P*{r) 



1/2 



r 



1/2 



(8) 

where we assumed the same mass density — and 
the same chemical composition over the entire scattering vol- 
ume 14- This is a reasonable assumption for a degenerate 



electron core. Near a SMBH, where orbital motion around a 
single mass dominates, the test particle (WIMP or star) ve- 
locities are Keplerian u» — v; in this case r/ = 3/2, al- 
though the exact value does not significantly change the re- 
sult. The value of the spin-independent WIMP-nucleon scat- 
tering cross-section ctq is limited by direct detection experi- 
ments, i.e. less than 10""'^ cm^ (CDMS 2006.) . If the star is 
composed of nuclei with atomic number An, the cross section 
increases by a coherent factor of A* . 

If a WD is heavy (M > Mq) and/or An > 1, al- 
most all WIMPs crossing the star will be captured. In this 
case, the WIMP capture rate is determined by the geometri- 
cal limit ttRI rather than the total interaction cross section 
(ToAnM^/Mn- We thus use a modified interaction cross sec- 
tion a'r, defined as 



Mn 



A 4 -^ -^^ 2 



(9) 



Figure[T]shows the capture rate by Oxygen WDs (An — 16) 
vs. WIMP velocity dispersion for several masses of WDs, 
assuming Keplerian orbits around the SMBH, — 100 
GeV, ao = 10-^3 cm\ ~ p™^'^ ^ mJ{{aav)T^i,) ^ 
10^° GeV cm^'^ which corresponds to the maximal cen- 
tral particle dark matter density allowed by the age of the 
SMBH Tbh ~ 10 Gyr and our selected value (agv) — 
3 X 10 -2^ cm3 s-i jBertone & Merrittl2005l;lGondolo & Sil5 
Il999l) . The left panel corresponds to WD effective tem- 
perature of T = 100, 000 K (without an envelope), where 
the masses and r adii used in the calculation are (Fig. 5 in 
iPanei et all 12000) : AU/Mq = 0.6, 0.8, 1.0, 1.2, 1.4, and 
i?*/i?Q = 0.02, 0.012, 0.0085, 0.006, 0.0045, coiTespond- 
ingly. The right panel corresponds to the zero-temperature 
approximation (Hamada & Salpeter 1961), where the mass- 
radius relation has been obtained by fitting the numerical re- 
sults for Carbon WDs with a function R^/Rq = 0.94 - 
O.67tan(1.49[Af*/M0 - 0.85]) in the interval AU/Mq = 
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Fig. 2. — Capture rate vs. distance from the central black hole for Oxygen (left) and Iron (light) WDs of T = 100, 000 K assuming Keplerian velocities for 
WIMPs and the WD. The WIMP mass density is normalized as (2 pc) = IOOMq pc~^ . For the central spike we assume a power-law profile wit h the indices 
(top to bottom): 7/3, 3/2, 4/3, where the maximal possible WIMP mass density p™'''' is given in the text. The mass-radius relation is taken from lPanei et all 
<25og) . 



0.15 — 1.4. The solid lines show the capture rate calculated 
using the modified interaction cross section (eq. ||9l), and the 
dashed lines are calculated for ctq = <to- In the case of Oxygen 
WDs the geometrical limit is reached for A/* ^ X.IMq; the 
larger mass WDs have smaller radii and therefore smaller cap- 
ture rates. For v < 10^ km s^^ our geometrical l imit calcula- 
tions agree well with the results of Bottino et al. (2002), their 
equation (26). For v > 10^ km s"^ the Bottino et al. (2002) 
formula, derived under the assumption that each WIMP cross- 
ing the star is captured, gives a systematically larger capture 
rate, up to a factor of 10 for 10^ km s^^. This can be treated 
as an upper limit, whereas our approximation to the geomet- 
rically limited case can be considered as a lower limit. 

It can be seen (Figure [U that cooler WDs have a capture 
rate (eqs. ||6l,||2l) larger than hot ones of the same mass be- 
cause of the larger escape velocity (eq. ||8]). The latter is 
the result of a smaller radius i?, and consequently stronger 
gravity. Thi s effect may b e explained in terms of the "focus- 
ing factor" ( lGouldllI987h or simply because WIMPs can be 
captured from the larger volume of the Maxwellian velocity 
phase space. A larger capture rate by a cooler WD will lead 
to accelerated heating until the WD radius increases due to 
increased temperature. A hotter WD will be less efficient for 
WIMP capture and cool down. This mechanism will thus lead 
to fast self-regulation of the WD temperature and capture rate. 

The capture rate for a different (An) composition WD can 
be estimated from scaling the Oxygen WD curves by a factor 
of (yl„/16)'^, e.g., the curves for a Carbon WD can be ob- 
tained from scaling the Oxygen WD curves down by a factor 
of (3/4)'^. WDs with heavier nuclei (up to iron) may exist 
jPanei et al.ll2000h : in this case, the capture rate is restricted 
mostly by the geometrical limit. 

Figure|2]shows the capture rate for Oxygen (left panel) and 
Iron (right panel) WDs vs. distance from the central black 
hole with Mbh = 3.7 x IO^Mq (iGhez et all [20051) : this 



includes effects of the radial dependence of the WIMP ve- 
locity dispersio n and the WD orbital (Keplerian) velocity. 
Follow ing Gnedin & Primacfl (i2004) and Bert one & Merrira 
dlOOl . the WIMP mass density is normaUzed as Pc) = 
IOOM0 pc~^. For the central spike we assume a power- 
law profile with indices (top to bottom): 7/3, 3/2, 4/3: these 
profiles are predicted for different scen arios of the black 
hole g rowth, adiabatic (|Ullio et ani2001l) . quasi- equilibrium 
( Gned in & PrimacS |2004|) . and instanteneous (lUlho et al.l 
i2001i) . correspondingly. Here we use the same estimate for 
^max ^Qj. pjgm-g [Y] f he mass-radius relation for Oxygen 
and Iron WDs of T = 100, 000 K is taken from lPanei etaP 
(.2000.) . 

The capture rate scales linearly with the WIMP density, 
so that the largest capture rate is reached with the adiabatic 
profile. The capture rate increases toward the SMBH until 
the maximal WIMP mass density p™^^ is reached; then the 
capture rate decreases due to increases in the WIMP velocity 
dispersion and the orbital velocity of the star. For the quasi- 
equilibrium profile, p™'^^ is reached only at ~10~^ pc, while 
the instanteneous profile is even flatter 

As can be seen from Figure |2] and a simple inspection of 
the capture rate formulae, in the geometrically limited case 
(eq. ||9l) the capture rate becomes essentially independent of 
the WIMP-nucleon scattering cross-section and degenerate 
core parameters. Observationally, the brightest WIMP burn- 
ers may be the geometrically limited ones. The main uncer- 
tainty in the geometrically limited capture rate is the dark mat- 
ter density: thus it may be possible to perform largely "model 
independent" measurements of the dark matter density pro- 
file by measuring the luminosity of different WIMP burners 
orbiting within a particular dark matter spike. 

A smaller annihilation cross section {<Tav) < 3 x 10"^^ 
cm^ s~^ would allow for larger ambient WIMP densities near 
the SMBH. This would lead to a larger capture rate and con- 
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sequently larger burning rate at the innermost radii. The en- 
ergy release due to WIMP annihilation in the stellar core is 
Ly, 0.16C(m^/100 GcV) erg s^^ which is actually inde- 
pendent of the WIMP mass m^. 

3. DISCUSSION 

Where does the ene rgy released during t he W IMP annihi- 
lation go? Table 1 in iMoskalenko & Wail (120061) shows that 
the effective radius of the thermal distribution of WIMPs in 
the stellar core is much smaller than the radius of the star 
<^ R*, therefore, the products of WIMP annihilation can- 
not propagate to the stellar surface and are converted into 
thermal energy and neutrino emission. A WD, a star with- 
out its own energy supply consisting of Carbon and Oxygen, 
may emit up to 3 x 10'^'* erg s^^, i.e. ^10 times the 

luminosity of the sun, burning WIMPs only and this energy 
source will last forever! (Note that this estimate is based on 
our approximation to the geometrically limited case and larger 
luminosities are possible even for the given set of parame- 
ters.) At such a luminosity, the surface temperature of the 
WD would be close to ~ 140,000 K, assuming = 1.2Mq, 
i?, = O.OO6i?0. The maximum of the black body emission 
falls into the UV band making such stars strong thermal UV 
emitters concentrated in the inner ^0.01 pc. A smaller an- 
nihilation cross section (aav) < 3 x 10~^^ cm'^ s~^ and/or 
larger WIMP density normalization (2 pc) > lOOMopc^"^ 
would allow for a larger ambient WIMP density near the 
SMBH, thus increasing the capture and burning rates further. 

The energy transport in the interiors of WDs is do minated 
by de generate electrons and is very efficient (see iHansenl 
|2004 for a recent review); therefore, the large number of 
captured WIMPs and their annihilation in the core would 
not change the internal structure of WDs. A recently pub- 
lished catalog of spectroscopically con firmed WDs from the 
Sloan Digital Sky Survey (SDSS) (E isenstein et"an l2006l) 
contains several hot WDs with surface temperature in the 
range of 100,000 K, thus providing observational evidence 
that high temperature does not change the appearance of 
WDs. A bare WD with an effective temperature as high 
as 170 ,000-200,000 K has also been observed (H1504H-65, 
IWerner & Wolff 1999). 

The number of very hot WDs in the SDSS catalog is small, 
just a handful out of 9316. This means that observation of a 
concentration of very hot WDs at the Galactic Center would 
be extremely unlikely unless they are "dark matter burners." 
The spectra of confirmed hot WDs can serve as templates for 
spectroscopic analysis of WDs at the Galactic Center where 
only a limited part of the near-IR band can be used. An inde- 
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pendent determination of the M^/R^, ratio is possible using 
the gravitational redshift that has to be equivalent to a ra dial 
velocity of about 50 km s~^ (iGreenstein & Trimbldl 19671) . 

A bare WD with a highly eccentric orbit around the central 
black hole may exhibit variations in brightness correlated with 
the orbital phase (Figure|2]i. To have this working, the orbital 
period should exceed the equilibrium time scale Teq (eq. JSj). 
Carbon burning stars hav e Tea ^ 10 yr, a nd it is even shorter 
-0.5 yr in case of a WD (IMoskalenko & Wai 2006). If a WD 
appears in a high-density WIMP region, the WIMP density 
in its material would quickly reach equilibrium; thus, the sur- 
rounding WIMP density variation as the WD orbits would re- 
sult in variations of brightness. This makes bare WDs ideal 
objects to test the WIMP density in the environment in which 
they are orbiting. Since cx p^, a population of WDs, 
bare or with envelopes, located at different distances from the 
SMBH would exhibit a luminosity correlated with the radial 
WIMP density profile. Geometrically limited WIMP burners 
have the highest luminosities and therefore will be the eas- 
iest to observe. Their luminosity is largely independent of 
WIMP-nucleon scattering cross-section, WIMP pair annihila- 
tion cross-section, and degenerate core parameters. 

Advances in near-IR instrumentation have made possi- 
ble observations of stars in the inner parse c of the Galaxy 
dGenzel et all 120001; iGhez et al.ll2003L l2005h . The apparent 
K-band brightness of these stars is 14-17 mag. The observed 
absorption line widths imply high temperatures and lead to 
a "paradox of youth:" apparently young stars in the region 
whose current conditions seem to be inhospitable to star for- 
mation. One of the possibilities is that they are old stars mas- 
querading as youths. Assuming a central spike with index 7/3, 
the K-band brightness for Oxygen WDs with T - 100, 000 K 
and i?*/i?0 ~ 0.01 is about 22-23 mag not including exti nc- 
tion, which may be as large as 3.3 mag ( Rieke et al.lll989l) . It 
is, therefore, unlikely that the currently observed stars in the 
K-band are WDs burning WIMPs; however, stars with degen- 
erate electron cores plus envelopes cannot be ruled out. 
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